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Abstract 

We present a detailed Monte Carlo study of the high-temperature phase of the 
two-dimensional driven lattice gas at infinite driving field. We define a finite- volume 
correlation length, study its finite-size-scaling behavior and extrapolate it to the 
infinite-volume limit by using an iterative extrapolation method. The same method 
is applied to the susceptibility. We determine the transverse exponents 7 and v. 
They turn out to be in perfect agreement with the theoretical predictions. We also 
compute the transverse Binder parameter, which apparently vanishes at the critical 
point, and the distribution of the magnetization. These results confirm the Gaussian 
nature of the transverse excitations. 



1 Introduction 



At present, the statistical mechanics of systems in thermal equilibrium is quite well estab- 
lished. On the other hand, little is known in general for non-equilibrium systems. It seems 
therefore worthwhile to study simple models which are out of thermal equilibrium. One of 
such systems is the model introduced at the beginning of the eighties by Katz, Lebowitz, 
and Spohn [![, who studied the stationary state of a lattice gas under the action of an 
external drive. The model, hereafter called driven lattice gas (DLG), is a kinetic Ising 
model on a periodic domain with Kawasaki dynamics and biased jump rates. Although 
not in thermal equilibrium, the DLG has a time-independent stationary state and shows 
a finite-temperature transition, which is however different in nature from its equilibrium 
counterpart 0. 

Despite its simplicity, the DLG has not been solved exactly ||. Nonetheless, many 
results have been obtained by means of Monte Carlo (MC) simulations and by using 
field-theoretical methods. In particular, Refs. developed a continuum model which 
is believed to capture the basic features of the transition and which provides exact pre- 
dictions for the critical exponents ||. Several computer simulations studied the critical 
behavior of these systems in two and three dimensions |7HT0| . These simulations provided 
good support to the field-theoretical predictions, once it was understood that the highly 
anisotropic character of the transition required some kind of anisotropic finite-size scal- 
ing (FSS). Some attempts to formulate a phenomenological FSS theory can be found in 
Refs. §p. 

In spite of the extensive numerical work, there are no direct studies of the correlation 
length so far, essentially because it is not easy to define it. Indeed, in the high-temperature 
phase the model shows long-range correlations due to the violation of detailed balance |12|, 
[13| . Therefore, no correlation length can be defined from the large-distance behavior 
of the two-point correlation function. A different definition is therefore necessary: in 



Refs. |14| a parallel correlation length is defined. However, this definition suffers from 



many ambiguities (see the discussion in Ref. and gives results for the exponent v\\ 
which are not in agreement with the theory 0. Even more difficult appears the definition 
of a transverse correlation length because of the presence of negative correlations at large 
distances |7|,|15|. 



In this paper we will define a finite-volume transverse correlation length generalizing 
the definition of the second-moment correlation length that is used in equilibrium systems. 
Because of the conserved dynamics, such a generalization requires some care. Here, we 
will use the results of Ref. . 



Once a finite- volume correlation length is defined, we can use the method of Ref. JT7 



to determine infinite- volume quantities |I8[. Such a method is particularly efficient and it 
has been applied successfully to many different models [1^-25]|. Once the infinite- volume 
quantities have been computed, we can determine the critical exponents by analyzing the 
singularities in the limit T — > T c . 

The paper is organized as follows. In Sec. ^] we describe the model and define the 
observables measured in the MC simulation. Section |3| reviews the FSS theory and intro- 
duces the basic formulae that are used in the analysis of the simulations. Next in Sec. |], 
we consider the field-theoretical model of Refs. MM in a finite geometry and compute the 



2 



FSS functions for several observables. In Sec. |5| we describe the simulations and present 
the results which are then discussed in Sec. || We confirm the field-theoretical predictions 
with good precision both for infinite-volume quantities and for the finite-size behavior at 
the critical point. Our results also indicate that the transverse Binder parameter van- 
ishes at the critical point, in agreement with the idea that transverse modes are purely 
Gaussian. 



2 Definitions 
2.1 The model 

As a prelude to the DLG consider a simple lattice gas of interacting particles. We consider 
a finite square lattice A and N particles, each of them occupying a different lattice site. 
A configuration of the system is specified by the set of occupation numbers of each site 
n = {rii G {0, l}}i S A- We assume a nearest-neighbor attractive ("ferromagnetic" in spin 
language) Hamiltonian 

H A [n] = -J 22 n i n h (!) 

where J > (in the following we will set without loss of generality J = 1) and the 
sum runs over all lattice nearest neighbors. We consider a discrete-time Kawasaki-type 
dynamics [EE] , which preserves the total number of particles iV or, equivalently, the density 



PA |A 



(2) 



ieA 



where |A| is the total number of sites in A. At each step, we choose randomly a lattice 
link If rti = rij, nothing happens. Otherwise, we propose a swap, i.e. a particle 

jump, with probability w(AH/T), where 

AH = H A [n'} - H A [n] (3) 

is the difference in energy between the new (n') and the old (n) configuration. If the 
probability w(x) satisfies 

w(—x) = e x w(x), (4) 

then the dynamics is reversible, i.e. satisfies detailed balance. Under these conditions 
there is a unique equilibrium measure given by the Gibbs one 

e -j3H A {n] 

= T e -pu K [nV ( 5 ) 

Z—/{n'} 

where f3 = 1/T. In the thermodynamic limit, i.e. |A| — ► oo, the lattice gas exhibits a 
second-order phase transition for p A = 1/2 at the Onsager value (3 C = |ln (y/2 + 1). 

The DLG is a generalization of the lattice gas in which one introduces a uniform (in 
space and time) force field pointing along one of the axes of the lattice, i.e. E = ESt: It 
favors (respectively suppresses) the jumps of the particles in the positive (resp. negative) 
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x-direction. If A is bounded by rigid walls, then E is a conservative field and it can be 
accounted for by adding a potential term to H\[n]. Therefore, the system remains in 
thermal equilibrium: the net effect of E is simply to induce a concentration gradient in 
the equilibrium state. 

Here, we will consider instead periodic boundary conditions ||27|| . In this case, the 



field E does not have a global potential and the system reaches a stationary state which, 
however, is not a state in thermal equilibrium. 

In the DLG transition probabilities take into account the work done by the field during 
the particle jump from one site to one of its nearest neighbors. In this case one proposes 
a particle jump with probability w(ft AH + ftEi) , with £ = (1, 0, —1) for jumps (along, 
transverse, opposite) to x. 

For E ^ and p\ = 1/2, the system shows a continuous phase transition at an inverse 
temperature /3 C (E) which saturates, for E — > oo, at ft c (oo) ~ 0.71/3 C (0). For ft < ft c (E) 
particles are homogeneously distributed in space, while for ft > ft c (E) the gas separates 
in two regions, one almost full and the other almost empty, with interfaces parallel to E. 

2.2 Observables 

We consider a finite square lattice of size L» x L± with periodic boundary conditions. We 
define a "spin" variable sj = 2n- } — 1 and its Fourier transform 



_ ,'2im 27rm\ 

k n , m = I "T 5 ~ r > \') 



where the allowed momenta are 



with (n, m) G Zl« X Zl ± . 

We consider the model at half filling, i.e. for p\ = 1/2. Then 

J> = 0, i.e. 0(k OiO ) = O. (8) 

jeA 

In the ordered phase |0(k)| takes its maximum for k = k 0) i, and the expectation value on 
the steady state of its module 

m(L||,L ± ) = ^(|0(k O)1 )|} (9) 

is a good order parameter. 

In momentum space the static structure factor 

G(k;L||,L ± )^(|0(k)| 2 > (10) 
vanishes at k 0j o because of Eq. (B) and is maximal at k 0j i, so that it is natural to define 



the susceptibility as [g£ 



X (L h L ± ) = G(koX,L h L ± ). (11) 
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Another interesting observable is the transverse Binder's cumulant g(L\\,L±) defined as 

9(i| ' ii) = 2 "lMPP' (12) 

Next, we would like to define a correlation length. In infinite- volume equilibrium systems 
there are essentially two different ways of doing it. One can define the correlation length 
in terms of the large-distance behavior of the two-point function or by using the small- 
momenta behavior of the two-point function (second- moment correlation length). In 
the DLG the first method does not work. Indeed, in the high-temperature phase the 
two-point function (s x so) always decays algebraically with the distance even in the high- 
temperature phase. Moreover, it is not positive definite because of negative correlations 



in the transverse directions ||15|| . This peculiar behavior is due to the fact that in the 
infinite-volume limit (at fixed temperature) the static structure factor G(k; oo, oo) has a 
finite discontinuity at k = 0. 

To overcome the difficulties of the real-space strategy we will define the correlation 
length by using the two-point function for small momenta. We follow Ref. |JL6[| , where we 
discussed the possible definitions of correlation length in the absence of the zero mode, 
as it is the case here. 

We consider the structure factor in finite volume at zero longitudinal momenta 

G ± (q;L h L ± )= G((0, q); L h L ± ), (13) 

(note that the conservation law implies G±(0; Ln, L±) = 0) and introduce a finite- volume 
(transverse) correlation length 



\ \G±(qj]L h L 



G±(qu L\\, L±) . 

11 1 - 1 , (14) 



where q n = 2 sin (irn/L±) is the lattice momentum. 
Some comments are in order: 

(i) If we consider an equilibrium system or a steady state in which correlations decay 
exponentially, then we have for q — > that 

Gl\q; L h L ± ) = X (L h L ± )-\l + ^(L h L ± ) q 2 + 0(q\ L~ 2 )). (15) 

where x{L\\, L±) is the susceptibility. Thus, £f,-(oo, oo) is a good definition of corre- 
lation length which has an infinite- volume limit independent of i and j. 

(ii) Since G±(Q; L\\, L±) = 0, qi and qj must not vanish. Moreover, as discussed in 
Ref. ||TB[, the definition should be valid for all T in finite volume. Since the system 
orders in an even number of stripes, for i even G±(qf, Ln, Lj_) = is zero as T — > 0. 
Therefore, if our definition should capture the nature of the phase transition, we 
must require i and j to be odd. Although any choice of i,j is conceptually good, 
finite-size corrections increase with i,j, a phenomenon which should be expected 
since the critical modes correspond to q — > 0. Thus, we will choose = (1,3). 



5 



For a more thorough discussion of the definition (14), which is mainly motivated by the 



absence of the zero mode, we refer to Ref. [16| where its finite-size behavior is studied in 



the context of the large- N approximation of the iV-vector model. 

Another quantity which is considered in the analysis is the amplitude A 13 defined by 



-2 



A 13 (L h L ± ) = ^-. (16) 

3 Finite-size scaling 

In the neighborhood of a critical point the behavior of long-range observables is controlled 
by few quantities, corresponding in the renormalizat ion-group language to coordinates 
parametrizing the relevant directions. When the system is finite, its size can be effectively 
incorporated in the scaling analysis as another relevant operator. This means that an 
observable O, which diverges in the thermodynamic limit as 

Oo(/3)~^ 70 fort = l-|--0+ (17) 

behaves in an isotropic finite system of linear size L as 

CM/3) « t-^F h0 (t- u /L) « L^F 2 , (t~ u /L) « L^ v F 3t o(U/L), (18) 

where £oo(/3) is the infinite- volume correlation length. From Eq. flT8[ ) we can derive a 
general relation for the ratio of Ol at two different sizes L and aL: 



o aL (P) „ ffa((3) 



(19) 



where we have traded £00/ L for fa/L by inverting fa « LF 3 ^(^ 00 /L). 

The function Fq(z) is universal and is directly accessible numerically, e.g., by MC 
simulations. It depends on a and, for a two-dimensional finite lattice of dimension Lj_ x L» , 
on the ratio L±/L». The appropriate FSS limit is indeed achieved by taking the limit 
L±, L\\,fa — > 00 with L±/L\\ and fa/L fixed. 

If we define z = fa/L, then z varies between and z*, where z* is defined by 



00). (20) 



The value z* is directly related to the behavior of the finite-size correlation length at the 
critical point, since fa{/3 c ) ~ z*L. For ordinary phase transitions z* is finite. At the 



critical point it also holds 
Then 

and therefore 



O l {P c ) oc Uf°' v . (21) 
FoV) = = oPol^ (22) 

70 = logFojz^) 
v log a 
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The above-presented results are valid for an isotropic system. On the other hand, there 
are clear numerical evidences that the phase transition in the DLG model is strongly 
anisotropic. The continuum field theories introduced to describe the critical singularities 
predict a nontrivial anisotropy exponent A. For example, the scaling form of the critical 
two-point function is predicted to be 

G{k ± , jfc||) « fi- 2+r, G{fik ± , / +A fe||), (24) 

where 77 is the anomalous dimension of the density field. 

It is then natural to assume the existence of two correlation lengths £±, £y which diverge 
with different exponents u± and v\\ related by 

v\\ = (1 + A)i/j_. (25) 

These considerations call for an extension of the FSS arguments. A phenomenological 



approach to FSS for the DLG has been developed |14| , keeping into account the strong 
anisotropy observed in the transition (for d = 2 see Refs. ||,|9|,[Tl|, for d = 3 see Refs. [PHI). 
Following this approach, we assume that all observables have a finite FSS limit obtained 
by taking the longitudinal size Ln and the transverse size L±_ to infinity keeping constant: 

• the anisotropic aspect ratio S = S& = L^ 1+A ^ / L±; 

• the FSS parameter £± !00 (P)/L± (or equivalently its longitudinal counterpart). 

Then, Eq. flTq) still holds by using the correct FSS parameter. Using transverse quantities 
we rewrite 

O l (J3) « t-^F h0 (r^/L ± ) « Lf^F 2j o(r^/L ± ) « L' r J ° /vx F 3 ,o(Zj_, 0O /L ± ). (26) 
Analogously Eq. ( |TP[ ) is recast in the form 

O aL± ((3,S) _ /£ ± , ix (/3) 



L± (f3,S) L 



a,S\. (27) 



where we have shown explicitly the dependence on a and on the aspect ratio S. Equation 
(|27| ) will be the basis for our analysis of the phase transition in the DLG . For the 
transverse finite- volume correlation length we will use £13 as defined in Sec. [2.2| . Moreover, 
as much of our efforts are aimed to test accurately the theoretical predictions of Refs. 
we assume that A = 2 for d = 2 and in the sequel we always consider S = S2. 



4 Field-theory description of the DLG 

As we explained in Sec. [2.1| , the DLG is a lattice gas model. However, in a neighborhood 
of the critical point (critical region) we can limit ourselves to consider slowly-varying (in 
space and time) observables. At criticality the lattice spacing is negligible compared to 
the length and time scales at which long-range order is established so that it is possible to 
formulate a description of the system in terms of mesoscopic variables. In principle, the 
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dynamics of such variables can be obtained by coarse graining the microscopic system. 
However, given the difficulty of performing a rigorous coarse-graining procedure, one 
postulates a continuum field theory that preserves all symmetries of the microscopic lattice 
model. By universality the continuum model should have the same critical behavior of 
the microscopic one. In this framework a field theory has been proposed m3|,[28[] (see also 
Ref. ||15|| ) and analyzed, giving rise to exact predictions for the critical exponents in space 
dimension d with 2 < d < 5. 



The starting point of the theory is the dynamic functional [|30|-|32"|j which reads (ne- 
glecting terms which are irrelevant by power counting) ||.^| 



J[s,S] = J d d xdt A |s[A _1 <9t + A ± (A ± - r) - pA||]s + iu V||Ss 2 + sA ± s| 



(28) 



where s(x, t) is the local density field (coarse-grained version of S{), s is the response field, 
t is the effective distance from the critical point, p is a parameter, and Uq is the coupling- 
constant of the theory which takes into account the leading effects of the microscopic force 
field. A careful analysis shows that only the parameter p needs to be renormalized, allows 
to determine the upper critical dimension, which turns out to be d c = 5, and gives exact 
predictions for the critical exponents 

For the structure factor the renormalizat ion-group analysis gives the scaling form 

G(k\\, All; t) = pr 2+r, G{k\\^- A , All//" 1 ; r^ 1/v ), (29) 

where, in d dimensions, 

V = 0, (30) 

v = I, (31) 



2 
1 



d), (32) 



and r oc T — T c . In two dimensions, for the transverse structure factor, this implies the 
simple scaling form 

G ± (k;r) = r- 1 f(k 2 /r). (33) 

Thus, in infinite volume we have 

& ~ t- v , X ~ r~\ (34) 

where v is given in Eq. fl3"T|) and 



7 = 1. (35) 

The function f(x) defined in Eq. (|33|) is trivial. Indeed, keeping into account causality 
and the form of the interaction vertex one can see that for fcn = there are no loop 
contributions to the two-point function (and also to the response function (s s x )). Thus, 
for all 2 < d < 5, Gj_(k, r) is simply given by the tree-level expression 

G^k;r) = ^-. (36) 
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Two observation are here in order. First, r is an analytic function of t such that r = for 
t = 0; thus, t = bt for t — > with b positive constant. Second, the function that appears 
in Eq. fl3"6|) refers to the coarse-grained fields, which, in the critical limit, differ by a finite 
renormalization from the lattice ones. Thus, for the lattice function we are interested in, 
in the scaling limit t — ^ 0, A; — ^ 0, with k 2 /t fixed, we have 

d±Mt{k]t) = ¥TVt> (37) 

where Z and b are positive constants. 

On the same footing, we can conclude that all correlation functions with vanishing 
longitudinal momenta behave as in a free theory. In particular, the Binder cumulant 
defined in Eq. ( |T2"D vanishes. It is also important to notice that Eq. ( |3"TD implies the 
exponential decay of 

G ±Mt (x ± ;t) = j d^ke 1 ^ G ±Mt (k;t), (38) 

which fully justifies our definition of transverse correlation length. 

In this paper, we will study the FSS behavior of the model. Here, we want to ana- 
lyze the corresponding continuum field theory in a finite geometry following the method 
applied in equilibrium spin systems (see, e.g., Ref. |33|, Chap. 36] and references therein). 
The idea is quite simple. Consider the system in a finite box with periodic boundary 
conditions. The finite geometry has the only effect of quantizing the momenta. Thus, 
the perturbative finite- volume correlation functions are obtained by replacing momentum 
integrals by lattice sums. Ultraviolet divergences are not affected by the presence of the 
box [34]] and thus one can use the infinite-volume renormalization constants. Once the 



renormalization is carried out, one obtains the geometry-dependent finite-size correlation 
functions. 

Following this idea, if we consider a finite box of size L» x L±, Eq. (|29|) becomes 

G(k\\,k±; r; L\\, L±) = yU~ 2 G(/c||/i _1 ~ A , k ± fi~ 1 ; T[i~ 2 ; L\\iT l ~ A , L ± ^l' 1 ), (39) 
which shows that at T = T c 

&(r e )~Lj_, X (T C ) ~ Lf~ L\. (40) 



Moreover, Eq. fl36|) holds in finite volume. Keeping again into account the relation between 
coarse-grained and lattice quantities, we obtain for the lattice correlation function in a 
finite volume in the continuum limit (i.e. in the FSS limit with k — > keeping k 2 /t fixed) 



Z(t;L h L ± ) 
k 2 + r(t;Ln,Lj 



G ±Mt (k;t;L h L x ) = u ^Jli % (41) 



where Z and r are analytic functions of their arguments. In the FSS limit, we expect 

Z(t;L,,,L ± ) = Z(tL 2 ± ,S A ), (42) 



T 



(t;L h L ± ) = L^_ 2 f(tLl,S A ). (43) 
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Using these expressions, in the FSS limit we find 
£i 3 (t;L h L ± ) 



((2n) 2 + f(tLl,S A )) 



-1/2 



.4 



13 



£l3(*;£||;-k-L) 

x(t]L h L ± ) 



Z(tL 2 ± ,S A ) 



(44) 
(45) 



valid for t — > 0, L»,L± — ► oo with 5*2 and tL^ fixed. Therefore, from the scaling of the 
correlation length and of the amplitude A13 we can derive the scaling functions Z and f. 

If we make the simplest approximations Z(tL 2 ± , 5a) = const and f(tLj_, 5a) = const x 
tL 2 ± , we obtain for the scaling functions defined in Eq. ( p?|) the approximate forms 



[1 - (1 - a" 2 ) (2.) V] 



-1/2 



[l-(l-a- 2 ) (2tt)V]- 1 . 



(46) 
(47) 



As we shall see, these expressions provide reasonably good approximations to our data, 
but do not describe the data exactly: The functions f and Z that are determined from 
the data are nontrivial. 

In the previous discussion we have neglected the possible presence of a dangerously 
irrelevant operator that becomes marginal at d — 2. Its presence may modify the scaling 
relations (p9|), fl33|), and (BIJ). Considering G±(k) we have 



G±(k; r; L«,Lj 



u 



[i 2 G±_ (kfj, ; r/i 2 ; L\\jx 



-l-A 



(4? 



where, a = (d — 2)/3 in d dimensions, and w is the irrelevant coupling. If the scaling func- 
tion vanishes for u — > 0, one obtains an anomalous scaling behavior. In two dimensions, 
since the operator is marginal (a = 0), we expect logarithmic corrections to the formulae 
previously computed. It is also possible that logarithmic terms modify Eq. (0). In the 
absence of any prediction, we will neglect these logarithmic violations. As it has been 
observed in previous numerical studies, if present, they are small 0. As we will discuss, 
this is confirmed by our numerical results. 

Finally, we want to notice that the Gaussian nature of the fields at zero longitudinal 
momenta allow us to write down the probability distribution of 0(k Ojn ). Indeed, the 
previous results give 



P(0(k o , n )) = ^exp 



L\\L 



± 



;-ko,n)G/ latt (k , n )<Kk , n ) 



(49) 



where N is a normalization factor. From this expression we can derive the probability 



distribution of M 2 = | 
Eqs. (|), (H), and (0)) 



(Ki)\ 2 /L 2 Ll 
We obtain 

P(M 2 )dM 2 



Of course, (M) = m, (M 2 ) = x/{L\\L±) (see 



- M2 /° 2 dM 2 



a' 



where 



a' 



UL 



-G 1 



,latt 



(2n/L_ 



UL 



(50) 
(51) 
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5 Numerical Simulation 



5.1 Setup 

We studied the phase transition of the DLG in two dimensions by MC simulations. Our 
aim was to check the validity of the FSS assumptions and eventually to use FSS to 
compute the critical exponents. The observables we took into account are those described 
in Sec. g^. 

We used the dynamics described in Sec. [O] with Metropolis rates, i.e. we set 

w(x) = min (l,e~ x ). (52) 

Simulations are performed at infinite driving field: therefore, forward (backward) jumps 
in the direction of the field are always accepted (rejected). 

The dynamics of the DLG is diffusive with dynamic critical exponent 2 = 4. Thus, 
it is important to have an efficient implementation of the dynamics in order to cope with 
the severe critical slowing down. 

We use the random-number generator ran_array described in Ref. |3f| in order to 
choose the links and to perform the Metropolis step and multi-spin coding techniques to 
evolve simultaneously many independent configurations. On 32-bit machines the state of 
a single site for 10 configurations is coded in a 32-bit word, using 3 bits per configuration. 
This peculiar allocation is useful when computing the number of occupied neighbors of a 
given site: Since this number is at most 4 in 2 dimensions we can use 32-bit arithmetics 
without overflow. The Metropolis step is coded in an entirely algebraic way avoiding 
branching instructions and taking full advantage of the pipelining capabilities of most 
microprocessors. On Pentium machines we achieved a speed of 0.01 /zsec for single spin 
flip, while our code performs a bit worse on Motorola i?10000 or Alpha processors but on 
these machines it is possible to simulate 21 systems using 64-bit integers. 

We considered lattices with an approximately constant aspect ratio S 2 ~ 0.200. We 
sampled very accurately those with sizes (L||,Lj_) corresponding to (21,14), (32,16), 
(46,18), (64,20), (84,22), (110,24), (168,28) for which S 2 lies between 0.197 and 0.202. 
We checked that a change of ±1 on L\\ has an effect on the observables that is negligible 
with respect to the statistical error. In the following we write L for L±. 

It is very important to be sure that the system has reached the steady-state distribution 
before sampling. Metastable configurations in which the chain could be trapped for times 
much longer than typical relaxation times in the steady state are a dangerous source of 
bias. In the DLG, configurations in which multiple stripes aligned with the external field 
are present are very long-lived and it is possible that they persist for times of the order 
of those of typical simulation runs, thus effectively inducing a spurious geometry on the 
system. To avoid the formation of stripes, we took care to initialize the larger systems by 
suitably rescaled thermalized configurations of smaller systems (where the stripes decay 
faster) at the same temperature and value of S 2 - 

We computed the autocorrelation time r x for the susceptibility x- Such an observable 
is expected to have a significant overlap with the slowest modes of the system, so that r x 
should give a good indication of the number of sweeps necessary to generate independent 
configurations. We found that, for the lowest temperature we considered (/3 = 0.311), r x « 
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L 


7 

412 


Cl3 









m 




3 








= 


0.2700 






14 
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0.23555(28) 


0.148869(50) 
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1.4448(15) 
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0.0786(13) 
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8.7094(53) 
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0.155093(51) 


0.22560(72) 
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1.5206(15) 


9.6826(74) 
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0.123347(51) 


0.16947(98) 


18 


1.5793(28) 


1.5796(18) 


10.4511(93) 


0.23873(41) 


0.100445(47) 


0.1273(11) 


20 


1.6280(35) 


1.6255(22) 


11.073(12) 


0.23862(49) 


0.082940(46) 


0.0923(13) 








= 


0.2800 






14 


1.5286(16) 


1.5172(11) 


9.4382(56) 


0.24389(27) 


0.161823(52) 


0.25090(68) 


16 


1.6159(21) 


1.6118(15) 


10.6473(81) 


0.24399(34) 


0.129594(53) 


0.19095(95) 


18 


1.6858(26) 


1.6821(19) 


11.641(10) 


0.24304(40) 


0.106178(49) 


0.1473(10) 


20 


1.7428(34) 


1.7408(22) 


12.421(13) 


0.24396(46) 


0.087965(48) 


0.1101(12) 








= 


0.2850 






14 


1.6086(16) 


1.5940(11) 


10.2434(60) 


0.24803(26) 


0.168977(54) 


0.27700(65) 


16 


1.7152(21) 


1.7048(15) 


11.7579(87) 


0.24719(33) 


0.136493(54) 


0.21965(88) 


18 


1.7998(27) 


1.7963(19) 


13.029(11) 


0.24765(38) 


0.112516(51) 


0.16941(98) 


20 


1.8654(34) 


1.8657(23) 


14.067(15) 


0.24745(45) 


0.093766(51) 


0.1334(12) 








= 


0.2900 






14 


1.6976(16) 


1.6769(12) 


11.1664(64) 


0.25183(26) 


0.176906(55) 


0.30542(61) 


16 


1.8243(22) 


1.8130(16) 


13.0604(95) 


0.25168(33) 


0.144208(57) 


0.24750(84) 


18 


1.9347(27) 


1.9256(19) 


14.729(12) 


0.25174(37) 


0.119924(54) 


0.20021(91) 


20 


2.0225(34) 


2.0137(24) 


16.134(17) 


0.25132(43) 


0.100638(55) 


0.1608(11) 


22 


2.0904(31) 


2.0862(21) 


17.310(15) 


0.25142(38) 


0.086542(40) 


0.1292(11) 


24 


2.1457(35) 


2.1447(23) 


18.304(17) 


0.25129(39) 


0.074305(35) 


0.1026(12) 


28 


2.2444(45) 


2.2389(29) 


19.900(22) 


0.25190(46) 


0.057904(33) 


0.0685(14) 








= 


0.2950 






14 


1.7941(16) 


1.7685(12) 


12.2087(68) 


0.25617(26) 


0.185533(57) 


0.33682(58) 


16 


1.9527(22) 


1.9338(16) 


14.602(10) 


0.25609(32) 


0.152967(58) 


0.28268(78) 


18 


2.0893(27) 


2.0740(20) 


16.797(14) 


0.25609(37) 


0.128498(56) 


0.23771(85) 


20 


2.2041(34) 


2.1941(25) 


18.754(19) 


0.25670(41) 


0.108808(60) 


0.1956(10) 


22 


2.2962(32) 


2.2878(22) 


20.434(18) 


0.25614(36) 


0.094271(43) 


0.1623(10) 


24 


2.3770(33) 


2.3731(23) 


21.937(19) 


0.25671(35) 


0.081524(36) 


0.1315(10) 


28 


2.5045(50) 


2.5028(33) 


24.333(29) 


0.25743(46) 


0.064138(40) 


0.0922(14) 








= 


0.2970 






14 


1.8368(16) 


1.8072(12) 


12.6676(70) 


0.25783(26) 


0.189260(58) 


0.35048(57) 


16 


2.0076(22) 


1.9852(16) 


15.298(11) 


0.25762(32) 


0.156795(60) 


0.29750(76) 


18 


2.1538(27) 


2.1369(20) 


17.709(15) 


0.25786(35) 


0.132092(58) 


0.25094(84) 


20 


2.2786(36) 


2.2698(26) 


19.949(21) 


0.25826(41) 


0.112359(62) 


0.2100(10) 


22 


2.3867(32) 


2.3778(23) 


21.940(19) 


0.25771(35) 


0.097791(45) 


0.1773(10) 


24 


2.4838(33) 


2.4787(24) 


23.735(20) 


0.25885(34) 


0.084924(39) 


0.1488(10) 


28 


2.6308(53) 


2.6194(36) 


26.590(34) 


0.25804(46) 


0.067096(45) 


0.1030(14) 



Table 1: Monte Carlo results. 



900, 1400, 2700 sweeps, for L = 20, 24, 28 respectively, where a MC sweep is conventionally 
defined as the number of moves equal to the volume of the lattice. For this reason, in 
order to have approximately independent configurations, we measured once every 1500 
sweeps. 

For each geometry, we performed long simulations at 15 different values of (3 in the 
disordered phase between 0.27 and 0.311, which correspond to temperatures ranging from 
1.41654 to 1.63217 in units of the critical temperature for the two-dimensional Ising model 
(3 C = log(l + v / 2)/2 ~ 0.440687. For each geometry and (3 we collected approximately 10 5 
measures. The raw data are reported in Tables [I] and |2|. 

The statistical variance of the observables is estimated by using the jackknife method [36]. 
To take into account the possible residual correlations of the samples, we used a blocking 
technique in the jackknife analysis. In the standard jackknife method the estimator for 
the variance is obtained discarding single data points. In the blocking technique, several 
variance estimators are considered, discarding blocks of data of increasing length and 
monitoring the estimated variance until it reaches a maximum. 
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L 


Sl2 


Sl3 


X 


A 13 


m 


a 




P = 0.3075 


14 


2.0812(17) 


2.0241(12) 


15.4022(79) 


0.26599(25) 


0.210342(60) 


0.42478(49) 


16 


2.3396(23) 


2.2947(17) 


19.719(13) 


0.26703(30) 


0.179677(65) 


0.38902(64) 


18 


2.5907(28) 


2.5477(22) 


24.273(19) 


0.26740(34) 


0.156237(65) 


0.35715(68) 


20 


2.8299(40) 


2.7938(30) 


29.116(31) 


0.26807(38) 


0.137182(79) 


0.32443(87) 


22 


3.0510(33) 


3.0174(25) 


33.872(28) 


0.26880(31) 


0.122863(55) 


0.30119(79) 


24 


3.2719(39) 


3.2349(29) 


38.911(36) 


0.26894(31) 


0.109908(55) 


0.27470(86) 


28 


3.6789(65) 


3.6355(49) 


48.776(74) 


0.27097(40) 


0.091858(75) 


0.2360(13) 




P = 0.3085 


14 


2.1082(14) 


2.0488(11) 


15.7157(69) 


0.26710(22) 


0.212671(52) 


0.43320(41) 


16 


2.3781(16) 


2.3270(12) 


20.2393(93) 


0.26754(22) 


0.182213(46) 


0.39818(45) 


18 


2.6418(18) 


2.5947(14) 


25.103(12) 


0.26819(22) 


0.159092(42) 


0.36907(47) 


20 


2.8958(21) 


2.8542(16) 


30.281(16) 


0.26903(22) 


0.140120(40) 


0.34012(51) 


22 


3.1328(36) 


3.0961(27) 


35.514(31) 


0.26992(33) 


0.125993(60) 


0.31663(82) 


24 


3.3744(39) 


3.3379(30) 


41.131(38) 


0.27089(31) 


0.113183(58) 


0.29252(82) 


28 


3.8316(69) 


3.7822(51) 


52.382(80) 


0.27308(40) 


0.095379(79) 


0.2569(13) 




P = 1.3000 


14 


1.8998(16) 


1.8659(12) 


13.3718(72) 


0.26038(25) 


0.194851(58) 


0.37052(54) 


16 


2.0932(22) 


2.0671(17) 


16.418(11) 


0.26026(31) 


0.162810(61) 


0.32164(73) 


18 


2.2586(27) 


2.2420(20) 


19.288(15) 


0.26061(35) 


0.138209(59) 


0.27788(78) 


20 


2.4132(36) 


2.3997(27) 


22.069(23) 


0.26093(40) 


0.118453(66) 


0.23766(98) 


22 


2.5425(33) 


2.5327(24) 


24.588(21) 


0.26088(34) 


0.103782(48) 


0.20613(97) 


24 


2.6602(34) 


2.6487(24) 


26.882(23) 


0.26098(32) 


0.090555(42) 


0.17492(97) 


28 


2.8636(55) 


2.8495(39) 


30.910(41) 


0.26269(44) 


0.072488(51) 


0.1300(14) 




P = 0.3025 


14 


1.9568(16) 


1.9155(12) 


14.0015(73) 


0.26205(25) 


0.199740(58) 


0.38791(52) 


16 


2.1682(22) 


2.1408(17) 


17.417(12) 


0.26314(31) 


0.168038(62) 


0.34245(70) 


18 


2.3616(27) 


2.3359(21) 


20.789(16) 


0.26247(34) 


0.143824(61) 


0.30275(75) 


20 


2.5389(37) 


2.5164(28) 


24.106(25) 


0.26269(40) 


0.124114(69) 


0.26489(95) 


22 


2.6910(33) 


2.6756(25) 


27.182(24) 


0.26336(34) 


0.109385(51) 


0.23325(95) 


24 


2.8355(35) 


2.8220(26) 


30.160(27) 


0.26404(32) 


0.096155(46) 


0.20347(94) 


28 


3.0806(57) 


3.0613(41) 


35.373(49) 


0.26494(42) 


0.077706(57) 


0.1565(14) 




P = 0.3050 


14 


2.0171(16) 


1.9683(12) 


14.6780(77) 


0.26395(25) 


0.204930(60) 


0.40679(51) 


16 


2.2498(23) 


2.2116(17) 


18.515(12) 


0.26418(31) 


0.173660(64) 


0.36476(68) 


18 


2.4701(28) 


2.4387(21) 


22.436(17) 


0.26508(34) 


0.149777(63) 


0.32770(72) 


20 


2.6757(37) 


2.6500(28) 


26.388(27) 


0.26613(38) 


0.130181(72) 


0.29203(89) 


22 


2.8641(34) 


2.8364(26) 


30.272(26) 


0.26577(33) 


0.115770(55) 


0.26642(88) 


24 


3.0339(36) 


3.0105(27) 


34.042(30) 


0.26624(32) 


0.102440(49) 


0.23553(89) 


28 


3.3461(65) 


3.3163(48) 


41.058(64) 


0.26786(44) 


0.083931(70) 


0.1886(14) 




P = 0.3100 


14 


2.1460(12) 


2.08097(94) 


16.1664(61) 


0.26787(19) 


0.215981(45) 


0.44418(35) 


16 


2.4343(12) 


2.37889(90) 


21.0314(69) 


0.26908(16) 


0.186047(34) 


0.41292(31) 


18 


2.7193(15) 


2.6677(11) 


26.3595(100) 


0.26998(17) 


0.163321(34) 


0.38646(37) 


20 


3.0030(17) 


2.9520(13) 


32.194(14) 


0.27068(18) 


0.144787(34) 


0.36129(40) 


22 


3.2738(37) 


3.2257(28) 


38.277(33) 


0.27184(33) 


0.131133(62) 


0.34229(79) 


24 


3.5552(41) 


3.5032(31) 


44.923(42) 


0.27319(31) 


0.118619(61) 


0.32141(81) 


28 


4.0953(72) 


4.0294(55) 


58.956(93) 


0.27539(40) 


0.101535(88) 


0.2944(13) 




P = 0.3105 


14 


2.1577(15) 


2.0938(11) 


16.3170(75) 


0.26867(23) 


0.217091(56) 


0.44857(42) 


16 


2.4548(17) 


2.3945(13) 


21.297(10) 


0.26921(22) 


0.187326(49) 


0.41820(45) 


18 


2.7459(20) 


2.6913(15) 


26.792(13) 


0.27035(23) 


0.164752(45) 


0.39204(48) 


20 


3.0396(23) 


2.9862(17) 


32.881(18) 


0.27121(23) 


0.146452(44) 


0.36941(52) 


22 


3.3226(37) 


3.2667(28) 


39.201(34) 


0.27222(33) 


0.132833(63) 


0.35071(79) 


24 


3.6092(41) 


3.5544(31) 


46.245(43) 


0.27320(31) 


0.120487(61) 


0.33219(80) 


28 


4.1856(73) 


4.1161(55) 


61.352(96) 


0.27615(39) 


0.103692(89) 


0.3059(12) 




P = 0.3110 


14 


2.1746(12) 


2.10673(94) 


16.4952(62) 


0.26906(18) 


0.218364(46) 


0.45238(35) 


16 


2.4750(17) 


2.4135(13) 


21.591(10) 


0.26980(23) 


0.188723(49) 


0.42338(45) 


18 


2.7757(20) 


2.7185(15) 


27.273(13) 


0.27097(23) 


0.166366(45) 


0.39995(47) 


20 


3.0736(23) 


3.0208(18) 


33.572(18) 


0.27182(23) 


0.148113(44) 


0.37754(51) 


22 


3.3632(38) 


3.3148(29) 


40.241(35) 


0.27305(32) 


0.134690(64) 


0.35948(77) 


24 


3.6772(41) 


3.6188(32) 


47.716(44) 


0.27446(30) 


0.122489(63) 


0.34152(78) 


28 


4.3024(76) 


4.2147(58) 


64.14(10) 


0.27697(39) 


0.106167(93) 


0.3206(12) 



Table 2: Monte Carlo results. 
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5.2 Finite-size scaling and extrapolation to the infinite-volume 
limit 



In this Section we want to test the validity of Eq. (|27j), by studying the FSS of the 
observables defined in Sec. |^. In order to compute F (z), we need to keep the ratio 
between the transverse sizes of two different lattices equal to a. The simplest thing to 
do would be to take a = 2 but this choice has the drawback of requiring large systems: 
indeed, if L± increases by a factor of two, L\\ increases by a factor of 2 1+A = 8. A more 
convenient choice consists in taking a noninteger a almost equal to 1, as far as it is allowed 
by the magnitude of the relative statistical errors. This strategy requires the knowledge 
of the values of the observables for noninteger values of L±. They are obtained by using 
an interpolation. In the analysis we used a = 1.25. 

At each (3 = (3 a , a = 1, . . . , n, our data consist of n a triples (Lf, £f , Of), i — 1, . . . , n a , 
of lattices and observed values. The data points are interpolated with smooth cubic 
splines £ a (L), O a (L). This interpolation provides estimates of £, a {L) and O a (L) for any 
(not necessarily integer) value of L between L min = min« Lf and £ max = maxj Lf. Then, 
we consider a set of m a = \ {n a + l)/2j values of L, 

= £ min + _ Tjj (« ~ L min) > ( 53 ) 

where j = 1, . . . , m a , and compute the corresponding triples 

3 3 U ' 3 \ L a 3 £*{Lf) O a (L]) J 

for j = 1, . . . ,m a . The triples ([54]) are the data we use in the computation of the FSS 
functions. Note that, since we need two independent data at different values of L to 
compute R, we generate a number of interpolated points that is one half of the total 
amount of data. This minimizes the correlations among the A a -. 

Errors are computed by using an auxiliary MC procedure that also estimates the 
covariance matrices (C|)jj = Cov (R^R^j) and (Cg)jj = Cov (Rq^Rqj). 

The points (zj,Rqj) are fitted to (z,F (z)) where F (z) is a polynomial in z of the 
form 

F a (z) = 1 + ao,i z 2 + a a ,2 z 3 + a a ,3 z 4 H (55) 

Such a parametrization, different from that used in Ref. |17[|, is partially motivated by 



the theoretical results of Sec. [4], indicating that, for small z, the FSS functions behave as 
1 + az 2 . If we neglect the correlations between z® and Rqj, then the x 2 minimization is 



a linear problem which is easily solvable [37 



In order to detect corrections to scaling, we perform the fit several times, including in 
each case only the results that satisfy Lu > L min C >, increasing step-by-step the parameter 

L m in,0- 

The number iVp arj e> of parameters in the fit as well as the value of -L m in,e> are chosen in 
such a way to optimize the values of the corresponding \ 2 - 111 particular, at fixed L m i n ,o, 
we take the value of N paXt o which minimizes x 2 /A^dof,e>) where A^ f,c» is the number of 
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^min ■ 


= 15.4 


-^min - 


= 16.8 


-^min 


= 18.2 


-^min 


= 21 






i? 2 (DOF) 




i?(DOF) 




iJ^(DOF) 




iJ^(DOF) 


0.27 


1.6703(69) 


1.2(2) 


1.6671(82) 


0.4(1) 


1.680(12) 


0.0(0) 






0.275 


1.8085(84) 


0.5(2) 


1.8042(95) 


0.3(1) 


1.820(13) 


0.0(0) 






0.28 


1.9827(98) 


0.9(2) 


1.979(11) 


1.0(1) 


1.996(17) 


0.0(0) 






0.285 


2.194(11) 


2.1(2) 


2.189(13) 


0.1(1) 


2.207(20) 


0.0(0) 






0.29 


2.477(11) 


3.1(5) 


2.470(13) 


3.2(4) 


2.495(17) 


1.3(3) 


2.502(25) 


0.7(2) 


0.295 


2.884(14) 


3.8(5) 


2.876(17) 


2.8(4) 


2.901(24) 


2.2(3) 


2.914(32) 


0.3(2) 


0.297 


3.090(15) 


6.1(5) 


3.082(18) 


5.8(4) 


3.101(26) 


6.2(3) 


3.108(37) 


6.4(2) 


0.3 


3.521(18) 


3.9(5) 


3.513(21) 


4.0(4) 


3.545(33) 


3.6(3) 


3.570(44) 


2.4(2) 


0.3025 


4.026(21) 


1.5(5) 


4.012(25) 


0.9(4) 


4.042(38) 


0.9(3) 


4.069(48) 


0.6(2) 


0.305 


4.793(25) 


10(5) 


4.790(30) 


4.6(4) 


4.816(41) 


1.9(3) 


4.854(58) 


0.4(2) 


0.3075 


6.269(33) 


3.4(5) 


6.257(38) 


3.4(4) 


6.256(52) 


2.4(3) 


6.332(72) 


0.05(2) 


0.3085 


7.391(38) 


7.7(5) 


7.390(42) 


4.7(4) 


7.357(58) 


2.9(3) 


7.464(85) 


0.7(2) 


0.31 


11.506(94) 


7.7(5) 


11.49(10) 


4.5(4) 


11.22(13) 


2.0(3) 


11.60(20) 


1.3(2) 


0.3105 


15.14(23) 


12.5(5) 


15.21(24) 


9.2(4) 


14.55(28) 


5.1(3) 


15.55(48) 


0.1(2) 



Table 3: Infinite- volume results for £ 13 , obtained from the extrapolation of finite-size data, 
for different values of L > L min . R 2 is the x 2 defined in Eq. flSTf), and DOF is the number 
of degrees of freedom. 



degrees of freedom of the fit (the number of data points minus iVp ari e>); then, we increase 
L m m,o until we attain x 2 I^Aoi,o ~ 1 an d the curve is compatible with the points with the 
largest value of L. 

In Fig. [I] we report the FSS function F^(z) for the correlation length £13. The data show 
quite a good scaling, although at a closer look one may observe some small systematic 
deviations. The results of the fits are illustrated in Fig. |2|. The figure is organized as 
follows. In the left frame we plot the residuals (the difference between the observed 
values and the fitting function) of the data used to perform the fit, i.e. the data points 
which fulfill the condition L > L m \ n ^; in the right frame the residual of all the data points 
with respect to the same fit. This is useful to detect possible systematic trends due to 
corrections to FSS that become smaller as the size of the system increases. The table 
under the frames reports the values of x 2 an d iVd f,£ for the fit as a function of the cut 
Lmm,£ an d of the number of parameters used N pSbT ^. 

We repeated the same analysis for £ 12 . As we discussed in Sec 2.2, this definition 
should not have a good behavior since it is not compatible with the symmetry of the 
low-temperature phase. This expectation is confirmed by the data, see Figs. || and |j: The 
results show strong corrections to scaling and the fitted function clearly underestimates 
the true Fg(z). 

Then, we considered the susceptibility as a function of £13/ L. The FSS plot is reported 
in Fig. |5| and the fit results in Fig. []. The corrections to scaling are larger than those for 
£13, but it is quite nice to observe that the results for the largest lattices coincide perfectly. 

As a further check we consider the ratio A13 defined in Eq. flTBp. Its FSS plot is 
reported in Fig. [7|. We observe here clear corrections to scaling that make impossible 
the determination of Fa{z). Note that, if 7/1/ = 2 as predicted by field theory, then the 
corrections must be nonmonotonic at least in the near-critical region. Indeed, if z* is the 
critical-point constant defined in Eq. (|20"D — as we shall see z* « 0.15 — , then Fa(z*) = 1 
if the exponents are the field-theoretical ones. On the contrary, our data seem to indicate 
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-^min - 


= 15.4 




= 16.8 


-^min 


= 18.2 




= 21 





Xoo 


i?(DOF) 


Xoo 


fi^(DOF) 


Xoo 


R z (DOF) 


Xoo 


R z (DOF) 


0.27 


12.049(42) 


0.8(2) 


11.958(49) 


1.5(1) 


13.9(1.7) 


0.0(0) 






0.275 


13.878(51) 


3.2(2) 


13.759(60) 


2.2(1) 


16.0(2.0) 


0.0(0) 






0.28 


16.352(64) 


0.8(2) 


16.210(75) 


0.4(1) 


18.9(2.3) 


0.0(0) 






0.285 


19.714(82) 


5.1(2) 


19.536(96) 


0.0(1) 


22.8(2.8) 


0.0(0) 






0.29 


24.683(87) 


7.9(5) 


24.49(10) 


5.0(4) 


29.1(3.4) 


1.3(3) 


26.03(79) 


1.3(2) 


0.295 


32.77(13) 


9.3(5) 


32.44(15) 


5.6(4) 


36.4(4.6) 


3.0(3) 


34.0(1.2) 


0.8(2) 


0.297 


37.49(16) 


15.5(5) 


37.10(18) 


11.1(4) 


38.6(5.1) 


9.3(3) 


38.8(1.3) 


8.7(2) 


0.3 


48.01(21) 


5.5(5) 


47.51(26) 


5.8(4) 


51.8(6.6) 


4.0(3) 


49.4(1.6) 


3.7(2) 


0.3025 


62.06(28) 


3.6(5) 


61.47(33) 


3.6(4) 


70.0(8.7) 


1.6(3) 


64.3(2.1) 


1.5(2) 


0.305 


87.02(42) 


14.0(5) 


86.35(48) 


8.4(4) 


91(12) 


4.4(3) 


89.8(3.0) 


0.8(2) 


0.3075 


146.30(83) 


4.8(5) 


145.01(95) 


4.9(4) 


158(20) 


3.4(3) 


150.2(5.1) 


0.1(2) 


0.3085 


201.4(1.4) 


11.2(5) 


200.0(1.5) 


7.9(4) 


211(27) 


5.7(3) 


205.8(7.1) 


1.2(2) 


0.31 


474.5(7.3) 


9.4(5) 


472.6(7.7) 


5.4(4) 


504(65) 


4.2(3) 


478(19) 


1.6(2) 


0.3105 


821(24) 


11.8(5) 


813(24) 


7.8(4) 


829(110) 


6.4(3) 


823(37) 


0.4(2) 



Table 4: Infinite- volume results for x, obtained from the extrapolation of finite-size data, 
for different values of L > L min . R 2 is the x 2 defined in Eq. flSTf), and DOF is the number 
of degrees of freedom. 



just the opposite: Apparently Fa(z*) is increasing with L and is larger than 1. Thus, we 
expect our FSS curves to be somewhat unreliable near z ~ z*\ This will be confirmed by 
a direct analysis that will be performed in the next Section. 

Finally, we report the FSS plots for the Binder parameter, Figs. ^, |^, and for the 
magnetization, Fig. [H]. The first one is quite good. There are tiny deviations, but they 
are not systematic as it can be seen from Fig. |9|. The somewhat large x 2 found in the fits 
is probably an indication that the statistical errors are slightly underestimated. The FSS 
plot of m is instead completely unreliable. At the values of L we consider, there are very 
large corrections to scaling and it is impossible to determine the FSS function F m (z). 

Let us now describe the extrapolation of £ and x to the infinite- volume limit. 

With our preferred values of L min and N pai we obtain estimates of the FSS func- 
tions and then, for each triple (L , £o(/3), Xo(P)) = (L,£({3, L), x((3, L)) obtained from the 
simulations we apply iteratively the map 

aL n _i \ 

u-mmn-m/Ln-!) , (56) 

Xn-l(0)4(£n-l/£n-l) / 

until the asymptotic limit (oo, £oo, Xoo(P)) is reached. Thus, for each finite-volume MC 
result 0(j3,L), we obtain an extrapolated value O 0O (/3,L). 

To combine the estimates C^^L,) with the same value of (3 and have a consistent 
estimate C 00 (/3), we proceed as in Refs. fI7|, We generate fake data sets, which 




are obtained by replacing each MC result O with error AO by O + rAO, where r is 
a normally distributed variable with zero mean and unit variance. Then, for each fake 
data set, using the same L miTlj0 and N parj0 , we compute the FSS curves and perform the 
extrapolations. Thus, for each (3 and L that we consider in the simulation we obtain a 
sequence of extrapolated values 0™(P,L), k — 1, . . . , M — typically M ~ 10 3 — where (3 



16 



and L indicate that 0£°(/3, L) is obtained by extrapolating the (fake) finite- volume result 
Ok(f3, L). We use these results to compute the covariance of the estimates at fixed (3 



1 



M 



CM 



M 



(57) 



k=l 



Then, we consider the weighted average 



0{(3) 



(58) 



which minimizes the residual sum of squares 



r 2 ((3) = j^ic-'imAOooi^u) - dmtPoofaLj) - dm. 



(59) 



If the variables C 00 (/3, Lj), i = 1, . . . , n(/3), are Gaussian with covariance C(/3), then 
R 2 {(3) is distributed as a x 2 random variable of n(j3) — 1 degrees of freedom (DOF). Then, 
the deviations of R 2 {(3) from n(/3) — 1 are a measure of the consistency of the different 
extrapolations and 0{j3) is our best estimate of Ooo{(3). The estimated error is 



The parameters of our preferred extrapolations are L min = 15.4 and iVp ar = 3 for £ i3 and 
L m m = 18.21 and iVp ar = 5 for \- The corresponding FSS functions are given by 



F x (z) = 1- 113.951 z 2 + 3053. 15 z 3 - 38192.7 z 4 + 204378 z 5 - 394893 z 6 . (62) 

The results are reported in Tables § and f|. All extrapolations are consistent, i.e. R 2 {(3) ~ 
n(f3) — 1, except in a few isolated cases (for example the extrapolation of £13 at (3 = 0.3105 
for L min = 16.8). Note also that the direct extrapolation of A 13 — here we use the larger 
lattices L > 21 to compute F^(z) — is in good agreement with the results obtained from 
the extrapolated estimates of x an d £, in spite of the fact that we do not have a reliable 
determination of Fa (z). Clearly the errors due the neglected corrections to scaling are 
small compared to the extrapolation errors. 

We wish finally to check the prediction flSDD for the probability distribution of M 2 . 
This cannot be done straightforwardly since we did not save the values of 0(k O) i) for 
each single configuration, but rather the sum of the values of 0(k o l ) for iV = 10 indepen- 
dent configurations — the ten configurations that were simulated together by our multispin 
program — . This means that we cannot determine from our data the distribution of M 2 , 
but rather that of 




(60) 




1 - 2.14297 z 2 - 171.098 z 3 + 535.11 z 4 , 



(61) 




(63) 



i=i 
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Pmin 


V 






X 2 /DOF 


DOF 


0.270 


0.556(15) 


0.552(22) 


0.31135(15) 


6.1 


11 


0.280 


0.542(13) 


0.577(21) 


0.31126(11) 


2.8 


9 


0.290 


0.531(11) 


0.600(21) 


0.311207(80) 


1.3 


7 


0.300 


0.518(15) 


0.629(32) 


0.311153(73) 


0.53 


4 


0.305 


0.522(56) 


0.62(14) 


0.31116(21) 


0.91 


2 


0.3075 


0.54(41) 


0.56(93) 


0.3112(13) 


1.0 


1 



Table 5: Results for the fit ([BSD for several values of /?, 



where 

M 2 L l{ L ± , , 

* = — , 64 

X 

where the sum is over N independent systems. Since the variables <7j are squares of nor- 
malized complex Gaussian random variables (see Eq. fl5"0|)), S is distributed as a reduced 
X 2 variable with 2N degrees of freedom (the two is due to the fact that the random vari- 
ables are complex). In Figure |TT]we report the distribution of S for some values of (3 and 
L, together with the theoretical prediction, that does not include any free parameter. As 
we can see the collapse is fairly good especially for the largest lattice near the critical 
temperature. 



5.3 Determination of the critical indices 

In order to determine the critical exponents we performed two different analyses. First, 
we considered the infinite-volume estimates of £oc(ft) an d Xooift) an d performed a fit to 
a power law in the reduced temperature t = 1 — /3//3 c . Then, we considered the FSS 
behavior determining the exponents from Eqs. (|2T|) and (f23|). 

Let us begin with the infinite- volume correlation length £ 13 (in the following we simply 
write it as £). We performed several fits of the form 

U(P)=B^t-\ (65) 

using in each case only the estimates of Ccx>(P) with (3 > m ™ - We find the results reported 
in Table [|. One clearly observes the presence of corrections to scaling: The estimate of 
v decreases as f3 m \ n increases until the fit looses predictivity. We may try to take the 
corrections into account, by performing a fit of the form 

U(P)=B^(1 + C^t 9 ). (66) 

Since no theoretical estimate is available for the exponent 9 and our data are not precise 
enough to keep 9 as a free parameter, we have performed several fits, using "reasonable" 
values of 9. 
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If we assume that the leading correction is the analytic one, i.e. 9 = 1, the fit improves 
substantially. The smallest value of x 2 /DOF is attained for the lowest possible value of 
Aniiu i-e. /3 m m = 0.27. Then, we obtain 

v = 0.4970(93), 

P c = 0.311109(44), 

Bi = 0.693(24), 

Q = -0.90(13), (67) 

with (x 2 /DOF) = 0.366 and DOF = 10. 

If we consider 9 = 0.5, we obtain again the best fit for /3 min = 0.27, yielding 

v = 0.459(16), 

P e = 0.311058(58) 

B ( = 0.866(63) 

G i = -0.658(92) (68) 

with (x 2 / D0F ) = °- 526 and DOF = 10. Note that the (x 2 / D0F ) is here slightly larger 
that in the case 9 = 1, providing some (very weak) evidence in favor of an exponent 9 = 1. 

These results indicate clearly that v < 1/2, with v decreasing with decreasing 9. If 
9 = 1, the result is perfectly compatible with the theoretical prediction fl3"4"|). We also 
tried to see the effect of an additional analytic correction. If we consider 

U% = B$r v (1 + C s t + D^t 2 ) , (69) 

we obtain for /5 min = 0.27 

v = 0.501(24), 

p c = 0.311118(73), 

B ( = 0.681(72), 

Q = -0.75(86), 

Dz = -0.6(3.8), (70) 

with (% 2 /DOF) = 0.399 and DOF = 9. As it can be seen from the errors on and D^, 
the addition of another fitting parameter is not justified. Our data are not precise enough 
to allow the determination of two correction-to-scaling terms. 

Finally, assuming the theoretical prediction v = 1/2, we compute our best estimate 
of the leading correction-to-scaling exponent. We redo the fit (|66D , fixing v = 1/2 and 
keeping 9 as a free parameter. Again, the best fit is obtained by using all data, i.e. with 
(3 min = 0.27. We obtain 9 = 1.09(22), p c = 0.311114(28) with (x 2 /DOF) = 0.356 and 
DOF = 10. Since an analytic correction should be present, and therefore 9 < 1, this 
result confirms the analytic nature of the first correction term. 

A further refinement can be obtained if we assume the theoretical prediction v = 1/2 
and 9 = 1. This provides our best estimate of the critical temperature. If we perform the 
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fit (|66|) fixing v and 9, we obtain 

p c = 0.311122(20), 

Bt = 0.6854(30), 

Q = -0.860(46), (71) 

with (x 2 /DOF) = 0.35 and DOF = 11. 

We can repeat the same analysis for the susceptibility \- We perform a fit of the form 



Xoo(P) = B x t~\ (72) 



and obtain for (3 min = 0.27 



7 = 0.995(34), 
(3 C = 0.31118(11), 

B x = 1.90(19), (73) 

with (x 2 /DOF) = 0.076 and DOF = 11. The x 2 is very small, indicating that probably 
the errors on Xoo(P) are overestimated. The final result for 7 is perfectly compatible 
with the theoretical estimate (|4"?D . Also, the estimate of /3 C is in good agreement with the 
estimates of (3 C obtained from the analysis of the correlation length. Although the goodness 
of the fit does not require it, we also performed a fit including an analytic correction. 
For (3 min = 0.27, we obtain 7 = 0.98(10), p c = 0.31116(20), (x 2 /DOF) = 0.0825 and 
DOF = 10. Clearly nothing changes. If we fix 7 = 1, we obtain a more precise estimate 
of p c . For j3 min = 0.27 we have 

P c = 0.311198(51), 

B x = 1.877(43), (74) 

with (x 2 /DOF) = 0.070 and DOF = 12. This estimate of f3 c is in good agreement with 
the result Q7Ip for the correlation length. 
We have also performed a fit of the form 

logx(^) = ailoge(/3) + 6 lj (75) 

that allows to determine a\ = j/u independently of the critical temperature. The results 
are reported in Table |6|. There are clear corrections to scaling, with the estimates of a\ 
increasing towards 2 as /3 m j n increases. These results are perfectly consistent with what 
we observed before. If no correction to scaling is included, the effective v varies between 
0.56 and 0.52 so that one would expect 7/V to vary between 1.8 and 1.9, which is exactly 
what we find here. 

In conclusion, we find good evidence that £ and x scale in the infinite-volume limit as 
predicted by field theory. If 7 = 1 and v = 1/2, we do not find evidence of nonanalytic 
corrections with 9 < 1. Assuming the theoretical predictions for the exponents we obtain 
for the critical temperature 

P c =0.31115(10), (76) 
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X 2 /DOF 


DOF 


0.270 


1.856(35) 


1.652(51) 


0.096 


12 


0.280 


1.867(42) 


1.631(66) 


0.10 


10 


0.290 


1.883(52) 


1.600(88) 


0.10 


8 


0.300 


1.899(51) 


1.57(10) 


0.041 


5 


0.305 


1.921(72) 


1.52(15) 


0.026 


3 


0.3075 


1.90(11) 


1.58(24) 


0.018 


2 



Table 6: Results for the fit (|75|) for several values of /?, 





w* 


70/1/ 


X 


1.543(2) 


1.944(5) 


A 13 


1.013(1) 


0.060(5) 


9 


0.902(4) 


-0.46(2) 


m 


0.804(7) 


-0.98(4) 



Table 7: Estimates of 70/^ for several different observables. 



where we have used the results ([71]) and ([74]) . The given error is a conservative estimate 
that should take into account both statistical and systematic uncertainties. The result 
for j3 c should be compared with the existing determinations: 

f 0.3108(11) Ref. [38]; , , 

Pc \ 0.3125(13) Ref. 0. 1 ' 

Our result (|76|) is in good agreement with both estimates, although much more precise. 

The critical exponents can also be computed from the FSS scaling curves. First, we 
determine z* by solving the equation a = F^(z*) and by using our determination of F^(z). 
We find z* = 0.1510(1). Then, for each observable, we compute w* = Fq{z*) and 70/V, 
see Eq. Q23p. The results are reported in Table [7]. For the magnetization m, cf. Eq. (^j), 
we extrapolated the data with L > 21, although they are probably not asymptotic. Thus, 
the correct estimate of /3/u could slightly decrease. 

The estimate of 7/1/ is not far from the expected value 7/1/ = 2. The result is slightly 
lower, which may indicate that the curve varies more strongly than predicted by our naive 
extrapolation of the FSS curve. The reason of the difference clearly appears if we consider 
v4i3 whose scaling dimension is 2 — 7/1/. Here, we find a slightly positive value related 
to the already observed fact that F^(z*) 7^ 1. Clearly for z«z* there are corrections to 
scaling that we have not taken into account. A better analysis will be presented below. 

Although our estimate of the FSS curve for m is quite bad, the result for (3/v is in 
very good agreement with the theoretical estimate (3/u = 1. Finally, let us consider the 
Binder parameter that apparently vanishes at the critical point as L~ 0A6 . The systematic 
decrease of the renormalized coupling can be seen more clearly in Fig. [1^. The fact that 
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7/1/ 


Z x 


x/dof 


DOF 


14 


1.975(24) 


0.0902(63) 


62 


5 


16 


1.955(24) 


0.0958(67) 


23 


4 


18 


1.939(27) 


0.1004(77) 


10 


3 


20 


1.925(31) 


0.105(10) 


3.207 


2 


22 


1.94(17) 


0.101(54) 


3.204 


1 



Table 8: Results for the fit (|8"0D for several values of L 



g scales with a negative exponent can also be obtained with a scaling plot, by looking for 
q so that 

gL q « F(£/L). (78) 
The optimal g is g ~ 0.5, in agreement with the result of Table [7|. The corresponding plot 



is given in Fig. 13 



It is also possible to compute the critical exponents by using the raw data near the 
critical point. The idea is to use Eq. (^1|), keeping into account the fact that our data 
are not exactly at (3 C . In this case we should add a correction term and use 

0{/3, L) = U' v (1 + a(P - (3 c )L l/u + ■■■), (79) 



where we have neglected corrections to scaling. Equation flT9|) is valid in the limit in 
which ((3 — (3c)L l l u is a small correction. 

We have considered our data at (3 — 0.311 that are quite near to the critical point and 
we have performed two sets of fits. In the first fit, we have neglected the correction term 
in Eq. ([T£]), by considering 

t(P,L) = Z 6 IS, X (P,L) = Z X L^. (80) 

Using our data for L = 14, 16, 18, 20, 22, 24, 28, we have performed a series of fits, including 
each time only the data with L > L mm . For £ the results are stable with L mm . The smallest 
X 2 is obtained for L min = 18. Correspondingly, 

p = 0.9923(85), Z ( = 0.1544(40), (81) 

with x 2 /DOF = 1.510 and DOF = 3. This is in perfect agreement with the expected 
value p = 1 and confirms the correct behavior of the correlation length. The analysis 
of the susceptibility gives the results reported in Table |8|. They show a clear downward 
trend and seem to indicate that 7/1/ < 2. This conclusion is however incorrect. The trend 
is due to the fact that we have not taken into account the fact that we are not exactly at 
the critical point. The correction term in Eq. (|79|) plays an important role, as we now 
show. We repeat the fit by using, in accordance with Eq. (|79|), 

X ([3,L) = Z X L^(1 + X X L 2 ), (82) 
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l/v 


X 


X x 


XVDOF 


DOF 


14 


2.060(43) 


0.0736(77) 


-0.000120(56) 


9.3 


4 


16 


2.031(80) 


0.079(16) 


-0.000091(90) 


7.5 


3 


18 


2.00(21) 


0.086(48) 


-0.00006(21) 


8.7 


2 


20 


1.9(1.1) 


0.11(34) 


0.0000(11) 


5.78 


1 



Table 9: Results for the fit (|82]) for several values of L mm . Here (3 = 0.311. 





l/v 




X x 


XVDOF 


DOF 


14 


2.061(50) 


0.0734(82) 


-0.000169(55) 


10 


4 


16 


2.037(85) 


0.078(17) 


-0.000147(90) 


9.5 


3 


18 


2.021(58) 


0.082(12) 


-0.000133(53) 


13 


2 


20 


1.9(1.2) 


0.11(37) 


0.0000(11) 


7.9 


1 



Table 10: Results for the fit (|82]) for several values of L min . Here (3 = 0.3105. 



where we have set v = 1/2 in the correction term. The results are reported in Table ||. 
They are now in good agreement with the theoretical prediction 7/z/ = 2. The same 
analysis has been performed for the correlation length. Using 

£(P,L) = Z^(l + X ( L 2 ), (83) 

for L min = 20 we obtain p = 1.000(94), = 0.1524(37), X € = 0.00(96) • 10" 4 with 
X 2 /(DOF) = 1.53, DOF = 2. In this case, the addition of the correction does not change 
the result. 

As a further check we have repeated the fit ( jB2|) by using the data at (3 — 0.3105. The 
results are reported in Table [H]. The estimates of 7/z/ are in good agreement with the 
expected 7/z/ = 2. Also, the estimates of X x are in rough agreement with the prediction 
X X (J3) oc {p - p c ), which implies X x (0.3105) w 4A X (0.3110). A combined fit of the data 
at p = 0.311 and (3 = 0.3105 to 

X ((3,L) = L^(1 + Y X ((3-(3 C )L 2 ), (84) 

gives, by using our estimate (|76|), j/is — 2.032(80) (resp. 2.00(22)) for L min = 16 (resp. 
18), with x 2 /DOF comparable to those of the previous fits. The agreement is clearly quite 
good. 

In principle, one should also consider corrections to scaling proportional to L~ w , where 
we expect uj > 2 as a consequence of our findings 9 > 1. However, our data are not precise 
enough to allow for the determination of such a term. 
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6 Conclusions 



In this paper we have performed a thorough check of the theoretical predictions for the 
DLG. Assuming A = 2, we checked that the transverse susceptibility x an d the transverse 
correlation length £13 have the correct behavior both in infinite volume for /3 — ^ /? c and 
at the critical point for L — > 00. In other words, we find very good agreement with the 
field-theoretical predictions ( [54]) and (f40|). Concerning the exponent /3, the analysis of 
the magnetization at the critical point gives, although with limited precision, (3 ~ 1/2, 
in agreement with the field-theoretical results. It is important to notice that in all these 
analyses we have not found any evidence for the presence of logarithmic corrections. As 
it has observed in previous studies ||, if they are really there, they are quite small. 

Our result for the Binder parameter g does not agree with that of Ref. || where 
it was found g 7^ at criticality, but confirms the results of Wang || that could not 
find a satisfactory collapse for the Binder parameter. Our result g = is compatible 
with the idea that in the scaling limit transverse correlations (both in infinite volume 
and in the finite-size scaling regime) are Gaussian, so that g — at criticality. Such 
a conclusion is also supported by the analysis of the distribution function of the finite- 
volume magnetization that is perfectly compatible with a purely Gaussian behavior. Note 
that this is not in contrast with the fact that g 7^ in the low-temperature phase. Indeed, 
it is possible that g behaves as the magnetization, i.e. g ~ (T c — T)£ for T — > T~, with 
C/ v = <? ~ 0-5- The observed exponent q is however difficult to explain. We mention that 
it is also possible that g decreases as a power of log L because of the marginal operator, but 
that, in our range of values of L, the complicated logarithmic dependence is mimicked 
by a single power. Note that, if gooiPc) = 0, the Binder parameter cannot be used to 
compute f3 c : The crossing method does not work. 

There is still a question that we don't fully understand. While the critical behavior of 
renormalized observables is fully in agreement with the predictions of the corresponding 
effective field theory, the FSS functions that we obtain for the microscopic model from 
our MC simulations are related to those of the coarse-grained effective field theory by the 
nontrivial scaling functions Z and f (see Eqs. ( fl~2"l) and (fE|)). In particular, the relation 
between the reduced temperature t of the microscopic model and the "temperature" r of 
the coarse-grained one seems to be not analytically predictable. Is this a consequence of 
the neglected marginal operator? An analysis of the three-dimensional case (where the 
operator is irrelevant), already in progress, should help us to gain a deeper understanding. 
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Figure 1: FSS plot for £ 13 versus £i 3 . The symbols correspond to the following lattice 
sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines connect 
data at the same temperatures. The dashed line is the approximation (0), the continuous 
line is the result of the fit with L mm = 15.4. 
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Figure 2: Fit of the FSS function of £13 versus £13. The symbols correspond to the 
following lattice sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. 
Dotted lines connect data at the same temperature. In the left frame only the data used 
in the fit are plotted, while in the right frame all the data are showed. The solid lines 
give the estimated error on the fitted function (one standard deviation). The graphs 
correspond to L min = 15.4. 
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0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15 0.16 

Figure 3: FSS plot for £ 12 versus £i 2 . The symbols correspond to the following lattice 
sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines connect 
data at the same temperatures. The dashed line is the approximation (fl6| ), the continuous 
line is the result of the fit with L min = 21. 
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Figure 4: Fit for the FSS function of £12 versus £12. The symbols correspond to the 
following lattice sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. 
Dotted lines connect data at the same temperature. In the left frame only the data used 
in the fit are plotted, while in the right frame all the data are showed. The solid lines 
give the estimated error on the fitted function (one standard deviation). The graphs 
correspond to L min = 21. 
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Figure 5: FSS plot for x versus £i 3 . The symbols correspond to the following lattice sizes: 
(o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines connect data 



at the same temperatures. The dashed line is the approximation (|47|), the continuous line 
is the result of the fit with L min = 18.2. 
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Figure 6: Fit of the FSS function of x versus £13. The symbols correspond to the following 
lattice sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines 
connect data at the same temperature. In the left frame only the data used in the fit 
are plotted, while in the right frame all the data are showed. The solid lines give the 
estimated error on the fitted function (one standard deviation). The graphs correspond 



33 




Figure 7: FSS plot for A 13 versus £13. The symbols correspond to the following lattice 
sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines connect 
data at the same temperatures. The continuous line is the result of the fit with L min = 21. 
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C(L)/L 

Figure 8: FSS plot for g versus £13. The symbols correspond to the following lattice sizes: 
(o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines connect data 
at the same temperatures. The continuous line is the result of the fit with L min = 21. 
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Figure 9: Fit of the FSS function of g versus £13. The symbols correspond to the following 
lattice sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines 
connect data at the same temperature. In the left frame only the data used in the fit 
are plotted, while in the right frame all the data are showed. The solid lines give the 
estimated error on the fitted function (one standard deviation). The graphs correspond 
to L min = 21. 
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Figure 10: FSS plot for m versus £13. The symbols correspond to the following lattice 
sizes: (o) L < 15, (□) 15 < L < 17, (A) 17 < L < 22, (0) L > 22. Dotted lines connect 
data at the same temperatures. The continuous line is the result of the fit with L min = 21. 
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Figure 11: Distribution function of the variable E. The symbols refer to different values 
of and L: o (/3 = 0.27, L = 16), □ (/? = 0.29, L = 28), A (f3 = 0.311, L = 16), 
(/3 = 0.311, L = 28). The solid line is the theoretical prediction as discussed in the text. 
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Figure 12: Binder cumulant: g versus £13/ X. The symbols correspond to the following 
lattice sizes: L= 14(B), 16(A), 18(T), 20(o), 22(D), 24(A), 28(0). 
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Figure 13: Rescaled Binder cumulant: L 5 g versus ^13/ L. The symbols correspond to the 
following lattice sizes: L= 14(B), 16(A), 18(y), 20(o), 22(D), 24(A), 28(0). 
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